Josephson current in unconventional superconductors through an Anderson impurity 
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Josephson current for a system consisting of an Anderson impurity weakly coupled to two un- 
conventional superconductors is studied and shown to be driven by a surface zero energy (mid-gap) 
bound-state. The repulsive Coulomb interaction in the dot can turn a vr junction into a -junction. 
This effect is more pronounced in p-wave superconductors while in high-temperature superconduc- 
tors with d^2_y2 symmetry it can exit for rather large artificial centers at which tunneling occurs 
within a finite region. 

PACS numbers;74.50.-hr,73.40.Gk,74.60.Jg 

Recently, the Josephson effect in unconventional su- 
perconductors have attracted a considerable attention 
10-^. Measurements of direct Josephson current yield 
valuable information on the symmetry of the order pa- 
rameter which is essential for understanding the mech- 
anisms of superconductivity in these complex materials. 
Phase interference experiments Q definitely suggest the 
presence of d-wave symmetry of the order parameter in 
high-Tc superconductors, while the recent discovery of 
superconductivity in Sr2RuOi Q implies the existence 
of a peculiar system for which the pair potential has a 
triplet (p-wave) symmetry 

p and d-wave symmetries of the order parameter have 
in common a property which reflects the variation of the 
pair potential across the Fermi surface. This results in 
a strong sensitivity to inhomogeneities which, in turn, 
influences the Josephson effect. In particular, an anoma- 
lous temperature dependence of a single Josephson junc- 
tion at low temperatures and an induced crossover 
from a usual (0-junctions ) to a 7r-junction on approach- 
ing the critical temperature were observed. 

Tunneling in a Josephson junction consisting of con- 
ventional (s-wave) superconductors and a dynamical im- 
purity (Anderson and Kondo ) was considered sometime 
ago |l0Hl^. In a recent work (l4j the tunneling current 
was calculated at zero temperature and was shown to be 
strongly dependent on the Coulomb interaction which, 
in some cases, may cause a sign change of the current. 
The experimental observations described above necessar- 
ily suggest that one must now go further and study the 
same device albeit with unconventional superconductors 
at finite temperatures. This missing part is addressed 
below. As will be demonstrated, the underlying physics 
is remarkably different. 

The main focus here is the influence of Coulomb inter- 
action on the low temperature behavior of the current in 
a 2D Josephson junction consisting of left (L) and right 
(R) superconductors with either p or d^2_y2 symmetry of 
the order parameter) weakly coupled to a quantum dot 
{via identical hoping matrix elements = Ir = t). The 



dot is represented by a finite U Anderson impurity whose 
energy — eo < lies below the Fermi energy. Usually, the 
inequalities U > —cq > Q are maintained so that the 
ground state of the disconnected {t — 0) dot is singly oc- 
cupied. We use the non-perturbative scheme suggested in 
Ref. [0 (extended for flnite temperature) and elucidate 
the low temperature behavior of the Josephson current 
and its dependence on U , t and the phase of the order 
parameter A. 

Since A is not isotropic, it is useful at this point to 
specify the underlying geometry. Each superconductor 
has the shape of half a plane defined as — oo < y < oo 
and a; < (x > 0) for the left (right) superconductor. 
The dot is located at the origin r = and tunneling is 
described by zero-range hoping between the impurity and 
the superconductors along the x axis. 

For d-wave superconductors we choose the nodal line 
of the pair potential on the Fermi surface to coincide 
with the tunneling direction, such that A = v^pp sin 2a 
where is the Fermi momentum. For spin-triplet super- 
conducting states the order parameter is an odd vector 
function of momentum and a 2 x 2 matrix in spin space. 
We chose to represent it by the time reversal symmetry 
breaking state Q which is off-diagonal in spin indices, 
that is, A = Aoexpia. This state is a probable candi- 
date for describing the recently discovered superconduc- 
tor Sr2RuOi §. 

In both cases, the pair potential has different values 
for electron like excitations which move at an angle a 
and hole like excitations propagating along the direction 
TT — a. This fact significantly affects the scattering pro- 
cess and causes the formation of a zero energy (mid-gap) 
bound state centered at the boundary. 

To obtain the Josephson current we compute the par- 
tition function Z = J Z?[?/;'0cc] exp(— S*) and the corre- 
sponding free energy of the system. The functional inte- 
gration is performed over Grassmann fields in the super- 
conductors and the quantum dot (see below for a precise 
definition). The Euclidean action can be written as a sum 
S = Sl + Sr -\- Sint + Sjj, where, with obvious notations, 
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(1) 

Sint^-T ^ (iiV';^i(0)r3C^ + i*cLr3V'c^j(0)), (2) 



ji=L,R 



Sl 



dT[cdrC + icT^c - U{ccf/2] 



(3) 



where e = eo + C//2. The summation is taken over odd 
Matsubara frequencies oj = {2n + 1)t:T, while ip^i = 
{ipuji1ipujii),c = (c-fC|) and the corresponding conjugate 
fields are Grassmann variables of the superconductors 
and the impurity respectively. 

The Bogolubov DeGennes hamiltonian ff,^^*^ ac- 
quires the form: 

ttBDG ^ ( e(p)-p|/2TO A(j5)exp(0,) 

where 4>i are the phases of superconducting condensates, 
and e{p) denotes the kinetic energy operator with disper- 
sion e{p). Note that we take = A/j but (pL ^ 4>R 

Since the pair potential is translation invariant in 
the y direction we can express the action in terms of 
fermion Grassmann variables V'ujfe (2;) and BDG Hamil- 
tonian H^^'~^{px,ky) where ky = ppsina. Moreover, 
due to the zero-range nature of the hopping, the fields at 
r ^ can be integrated out yielding effective actions in 
terms of the boundary fields "0^^5,(0) alone iQ. Explic- 
itly, 

Sl=TY, i^Luk, imLkjLuk, (0), (4) 

LJ ,ky 



with, 



GLky^^f [ dxdx'T^GLkAxx'), (5) 
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The Green func- 



where ko = ypp ,.y. 

tion Ghkyixx') satisfies the equation — T3[ia; + 
HE°'^{Px,ky)]GLky{xx') = S{x - x') with homoge- 
neous boundary conditions. It is constructed as a 
sum over the basic solutions of the equation —t^Iilu + 
HE^'^{px,ky)]uLky{x) — with coefficients depending 
on x' . The spinors u^kyix) describe electron like and 
hole like excitations moving in different pair potentials. 
For p-wave superconductors one has GLky — 2™I;''i(^) 
with 
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(6) 



It is markedly different from the corresponding quantity 
derived for s-wave superconductors (note in particular 



the occurrence of uj in the denominator) . Performing the 
integration over the boundary fields 'ipLujky{^),i'Luiky{^) 
in the partition function is now straightforward. The 
remaining integrations over the c, c fields is done by 
decoupling the Hubbard interaction using a Hubbard- 
Stratonovich transformation. Since the impurity oper- 
ators do not depend on ky, the auxiliary fields 7 depend 
only on w. We then obtain 



Z^C J];d7^exp(- 



2UT T 



(7) 



Here Q — —Tj^^j ln,[detR{LLj)] and C is a constant. The 
matrix R encodes the coupling between the supercon- 
ducting (half) planes connected and the impurity, and 
reads 

R{lu) = iT3 + 7^ - - ^^^^ < \t\^{rL{uj) + TRicu)) >, 

UJ 

(8) 

where N{Q) is the density of states at the Fermi level, 
and < O >= TT^^ \'^^i^dctO(a). The averaging proce- 
dure thus defined is related to a possible dispersion of 
the transmission matrix element t^y. For a point junc- 
tion with constant t, only the component of the Joseph- 
son current perpendicular to the interface is relevant. In 
this case, for even (d-wave symmetry) superconductors 
the mid-gap zero energy bound state does not contribute 
to the Josephson energy. Contrary, for p-wave pair po- 
tential this state defines the main contribution at low 
temperatures (which is the temperature domain of 
our interest here). Intuitively, occurrence of dispersive 
tunneling matrix elements t^j, correspond to deviation 
of the impurity from a point-like defect. Such a case 
can be realized e.g. by artificially - induced defects p6| . 
The spectroscopy of BiiSriGaCuiO% surfaces indicates 
that such defects appear to be more extended in STM 
imaging. In this case one can expect non-zero contribu- 
tion from the mid-gap level in d-waves superconductors 
as well. 

With this point in mind we now proceed and consider 
the p-wave case. The functional integral in (7) is approx- 
imated by the saddle point method. The appropriate 
optimum solutions 7cj should then minimize the free en- 
ergy 



F = -T 



(9) 



Here we use the notations: 



A{u;) = + co^il + a{u;)r - 7^ - b\u;)\A\^ cos\^), 

(10) 



2 



a{bj) ~ r 



(11) 



where F = 27r<^A^(0) is the bare impurity level width, 
b{uj) = T/lu , and (j) is the phase difference between the 
two superconductors. The self-consistent equation for 
is 



aj>0 



= 0. (12) 



Once the solutions are defined we can calculate the cur- 
rent J = (2e/h)dF/d(j) and the impurity occupancy 
n = dF/deo, 



J 



-(2e/?i)sin(0)r^ 



A\%^{u;)A{uj) 



w>0 



i-4r^ 



w>0 



A2(cc;)-t- 472^2(1 ^„(^))2 



(13) 



(14) 



We now analyze the main results of the present study. 
All the parameters having the dimension of energy 
(e, J7, r,T) are expressed in units of |A| and the current 
is given in units of \A\e/h. 
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In Fig. la, b, the current is displayed versus temperature 
in the low temperature region and for coupling strengths 
F ranging between 0.001 and 0.0882. At small Hubbard 
interaction (J7=2.1, Fig. la) the junction is in a V-state 
for which the current (within the present geometry of the 
tunneling direction) is negative. The current is strongly 
dependent on temperature, which is markedly distinct 
from the classical Ambegaokar-Baratoff formula. It is 
typical for superconducting systems for which the zero 
energy mid-gap bound state plays the major role. At 
higher values of U (Fig. lb) the pattern is inverted: The 
current for most coupling strengths is now positive and 
the 'tt' junction is transformed into a '0' junction with 
slightly weaker temperature dependence. A similar situ- 
ation takes place also in s-wave superconductors ||l2|,|l^ , 
where it is attributed to the single occupancy of the im- 
purity which then becomes a degenerate magnetic mo- 
ment. We have thus presented another example for this 
scenario, though the temperature dependence is quite dif- 
ferent. 
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FIG. 2. Dependence of the Josephson current on the phase 
difference for (a) small transparencies and intermediate tem- 
perature (r = 0.01) and (b) large transparencies and very low 
temperature (T = 0.001). Other parameters are eo = —2.0 
and U = 2.6. 
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FIG. 1. Dependence of the Josephson current on tempera- 
ture for U = 2.1 (a) and U = 3.0 (b) for different transparen- 
cies P. Note the sign reversal for strong U. Here eo = —2.0 
and (j) = 7r/12 are fixed. The units of all quantities are ex- 
plained in the text. 



Now let us discuss the behavior of the current J{(t>) as 
function of the phase difference as displayed in Figs. 2a, b 
for fixed U and T and for numerous transparencies F. 
Once again it is clear that as the transparency is in- 
creased, the sign of the current is reversed. Moreover, at 
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lower temperatures (Fig. 2b) when F/A > 0.4 it becomes 
irregular with jumps at certain values of </> at which the 
impurity is nearly incompressible. This effect can mani- 
fest itself when the flowing current I > Ic- In the frame 
of the RSJ model the averaged voltage V on the junction 
is related to the current / and the resistance R as, 

Jo 2n{I-J{^)) 

The irregular pattern of J((/>) results in a deviation 
of the V{I) characteristic from its classical expression 
/ = ^/J^ + [V/Rf. Such deviation is seen on Fig. 3 
where the difference is mainly exhibited at low voltage. 
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FIG. 3. V(I) characteristics of the p wave junction 
(upper curve) compared with the classical square root 
expression (lower curve). The relevant parameters are 
eo = -2.0, U = 3.0 and T = 0.66. 

To summarize, we have solved the problem of transport 
between two unconventional superconductors through 
an impurity and traced the dependence of the cur- 
rent on temperature. Coulomb interaction barrier trans- 
parency and phase of condensates. The essentially 
non-perturbative, self-consistent approach we have used 
yields a flnitc value for the current without adding re- 
laxation terms that should be included if perturbation 
approach is adopted. It is mandatory at this final stage 
to point out the peculiarities related to the physics of non 
s-wave superconductors. 

1. The contribution to the current in our case originates 
principally from the surface bound state which is related 
to the asymmetry of the pair potential, and has no ana- 
log for s-wave superconductors. The contribution of this 
state in low transparency junctions results in a large cur- 
rent, that is, \Jp/Js\ ~ -\/A/F. It is also marked by a 



stronger temperature dependence especially at low tem- 
peratures. 

2. For unconventional superconductors, the Josephson 
tunneling through an interacting quantum dot can serve 
as an indicator to distinguish odd parity supercondiictors 
from even parity ones. Superconductors with both p and 
d-wave symmetry of the order parameter have surface 
bound state which contributes to the current (mainly at 
low temperatures). Yet, as we indicated above, for even- 
symmetry (d-wave) superconductors, the current van- 
ishes in the limit of point-like impurity. In sharp distinc- 
tion, under the same conditions, the current is maximal 
for odd-symmetry (p-wave) superconductors. 
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